addpath ../Generator
addpath ../Generator/GrTheory
addpath ../PDCO
addpath ../CRNTSolvers


%load 'SimpleNetwork.mat'
n = size(Ak,1);
m = size(Y,1);
%Establish an initial value for the mu parameter
x0 = ones(n,1);
%call the solver
[iter,v_vecs,lmda_vecs,mass_infeas,mass_action_infeas,mu_vecs]=...
		FpIterationPDCO(Y,Ak,1e-6,500,x0);
%Calculating the concentrations of the speces
C = exp(lmda_vecs);
plot(C');
title('Concentration paths of species vs iteration of fixed point algorithm');
ylabel('Concentration of species');
xlabel('Iteration');
